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We investigate the hadron-quark phase transition inside neutron stars and obtain mass- 
radius relations for hybrid stars. The equation of state for the quark phase using the 
standard NJL model is too soft leading to an unstable star and suggesting a modifica- 
tion of the NJL model by introducing a momentum cutoff dependent on the chemical 
potential. However, even in this approach, the instability remains. In order to remedy 
the instability we suggest the introduction of a vector coupling in the NJL model, which 
makes the EoS stiffer, reducing the instability. We conclude that the possible existence 
of quark matter inside the stars require high densities, leading to very compact stars. 
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1. Nambu-Jona-Lasinio model 

A study of the stability of quark matter was carried on in RefP using the Nambu- 
Jona-Lasinio model in the SU(2) version with a repulsive vector coupling. In that 
work the pressure and the energy per particle, as a function of the baryonic den- 
sity, was analyzed. We also discussed the influence of the vector interaction in the 
equation of state (EoS) and studied quark stars composed of pure quark matter 

with two flavors. In the present work we consider the same NJL Lagrangian in the 
forrrPS 

C = - m)q + g s [{qqf + {&K?qf] ~ gviqi^q) 2 , (1) 

where q is a fermion field with Nf = 2 flavors and N c colors. Except for the bare 
mass m, the Lagrangian is chirally symmetric SU(2)l x SU(2)r. We included in- 
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teraction terms in the scalar-isoscalar, pseudoscalar-isovector and vector-isoscalar 
channels. The gs and gy are the scalar and vector couplings, respectively, and they 
are assumed to be constants with dimensions MeV~ 2 . 

Expanding qq and qj^q we can derive the mean field thermodynamic potential 
at temperature T and chemical potential //. We restrict ourselves to the Hartree 
approximation. The thermodynamic potential density, J7, depends on two parame- 
ters, namely the dynamical fermion mass, M, and the renormalized quark chemical 
potential, /ir, which are related to th e scalar, (^ip), and vector, (ip^ip), densities at 
the chemical potential /i throug 

M = m-2ga(H), ( 2 ) 

M = li-2gv(il>H)- (3) 

These are the NJL gap equations for the dynamical mass and the renormal- 
ized chemical potential. The vacuum contribution to the thermodynamic potential 
density is divergent and has to be regularized. This regularization is performed by 
introducing a cutoff, A. Once we know the thermodynamic potential density, other 
thermodynamic quantities such as the baryon number density, the energy density 
and the pressure, can be calculated in the standard way-M 



2. Hybrid Stars Phenomenology 

In the interior of astrophysical compact objects such as a neutron star, the density 
of matter can be several times the nuclear matter saturation density. Calculations 
based on microscopic equations of state, which include only nucleons as degrees 
of freedom, show that the central density of the most massive neutron stars is 
^7—10 times the saturation density of nuclear matter, pq. The precise point where 
the transition from nuclear matter to the quark matter occurs, is not known and 
only estimates have been put forward in the literature! We can assume a first order 
phase transition, as suggested by lattice calculations, and use the mentioned NJL 
model to describe the quark phase and estimate the phase transition within the 
compact object^. This approach has been suggested by several authors and there is 
a vast literature in this topic, including calculations wit h ot her relativistic models 
such as the MIT and the chromodieletric model (CDM^^-A Neutron stars having 
deconfined quarks in their interior are known as hybrid stars. 



2.1. Phase Transition 

According to the Gibbs criterium, two phases are in thermodynamical equilibrium 
when their baryonic chemical potentials, temperatures and pressures are equal, 
corresponding, respectively, to chemical, thermal and mechanical equilibrium^: 

MH = MQ = M; (4) 



Th = T q = T, 



(5) 
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and 

Ph(M,T) = P q ( M ,T)=P. (6) 

Here, the indices H and Q refer to the hadron and quark phases, respectively. For 
the hadron phase we use the equation of state known as GM1 (Glendenning and 
Moszkowskij^l and for the quark phase, we use the mentioned NJL model with 
vector coupling. The GM1 parameterization describes the hadronic matter in the 
interior of the star and the NJL model the quark matter star core. 

The EoS for quark matter obtained in RefP is used to solve the Tolman- 
Oppenheimer-Volkoff equations. These equations are integrated from the center 
of the star, where M(R = 0) = and e(R — 0) = e c , with e c the central energy 
density, up to r — R, where the pressure drops to zero. For each E oS we o btain a 
unique relation between the star mass and its central energy 

density^EKO. Fig.rjJ 

illustrates the behavior of the mass of the neutron star as a function of the radius 
and of the central density for different ratios gv/gs- We observe that the mass is 
characterized by a "cusp" which, according to Baldo et alW, is related to an in- 
stability. The plateau that appears in Fig. [1] (right) is a consequence of the Gibbs 
criterium. The width of the plateau is related to the jump in the density at the 
onset of the quark matter in the interior of the star. Baldo et al., based on the 
fact that the quark EoS using the NJL model is too soft to yield a stable star with 
quark matter in its interior, suggested a modified NJL model with a momentum 
cutoff that depends on the chemical potential. However, even in this approach, the 
instability remainJ^l. Therefore, we suggested the introduction of a vector coupling 
in the NJL model, making the EoS predicted by this model more stiff, as discussed 
in Refffl Hence, including the vector term in the Lagrangian density, obtaining the 
EoS for a quark phase and taking into account a hadron phase described by the 
GM1 approach, we can investigate a possible phase transition, finding, as Fig. [T] 
shows, that by increasing the vector coupling, the "cusp" that appears in the mass- 
radius graph, becomes more tenuous, and the plateau that appears in the mass 
versus central density plot decreases from gv 1 9s — to gv/gs = 0.25 and even dis- 
appears at gv/gs = 0.5. We do not show the effect of higher couplings, because, in 
that case, the phase transition starts at densities much higher than those considered 
here (p < 10po)- 



3. Conclusions 

For a vector coupling such that gv/gs = 0.5 we observe that the quark phase 
shows up right before the second maximum in the plot of the mass versus central 
density (the small central density range where d^/-^™) > 0). We know that 
hadronic matter is stable up to the first maximum. The fact that there is a second 
maximum, means that in a small range of central densities, p c /po ~ 9, the star is 
stable. The smallness of this region indicates that the possibility for stable hybrid 
stars is remote. We conclude that for the phase transition from nuclear to quark 
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matter to take place for larger values of gv/gs requires very high central densities, 
leading to more compact stars. Finally, from this study on hybrid stars, we conclude 
that the use of the NJL model with vector coupling restricts the threshold value for 
the vector repulsion. 
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